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SECTION  I 


INTRODUCTION 


SECTION  I 


Introduction 

The  object  of  this  technical  report  'is  to  inform  the  Office 
of  Naval  Research  on  the  present  status  of  their  Project 
NR  064  369,  Contract  No0  N7onr-32911 , "The  Dynamic  Properties 
of  Plastics  and  Rubber-like  Materials",  which  is  being  con- 
ducted by  the . Department  of  Mechanics  at  Illinois  Institute  of 
Technology,, 

Two  apparatuses  are  being  used  to  study  the  effects  of 
dynamic  loading  on  materials: 

The  first  uses  mechanical  and  optical  devices  to  determine 
directly  the  stress-strain  curves  of  plastic  and  rubber-like 
materials  subjected  to  impact  loads,  the  durations  of  which  are 
of  the  order  of  milliseconds. 

The  second  apparatus,  which  is  only  partially  built,  will 
use  mechanical  and  electronic  devices  to  determine  the  stress- 
strain  curves  of  elastic  materials  subjected  to  impact  loads, 
the  durations  of  which  will  be  of  the  order  of  microseconds. 

In  this  report  a general  description  of  the  apparatuses  is 
presented  and  the  theories  of  both  methods  are  outlined  briefly. 
The  mechanical  and  optical  parts  of  the  apparatus  for  testing 
materials  at  low  rates  of  compressive  straining  are  described, 
the  theory  of  the  experimental  method  is  outlined,  and  the 
influence  of  the  principal  causes  of  mechanical  and  optical 
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error  associated  with  the  method  are  evaluated,,  It  is  shown  how 
the  dynamic  stress-strain  curves  for  a given  specimen  under  im- 
pact loading  can  be  calculated  directly,,  when  the  hereditary 
characteristics  of  the  material  are  known . By  comparing  the 
experimental  and  calculated  curves  it  is  possible  then  to 
determine  the  accuracy  of  this  experimental  method  of  in- 
vestigating the  dynamic  properties  of  materials „ Some  results 
are  given  for  a natural  rubber  of  30-durometer  hardness. 

The  component  mechanical  devices  and  electronic  equipment 
for  testing  materials  at  higher  rates  of  compressive  straining 
are  described  and  the  theory  of  the  method  is  given. 

A mathematical  analysis  of  the  longitudinal  impact  of  a 
round- he aded,  short  bar  on  an  infinitely  long  bar  is  presented. 
This  study  follows  that  given  by  W,  A.  Prowse  [1],  C2],  [3],  [4], 
but  the  analysis  is  extended  so  as  to  give  final  results  in  an 
explicit  form  more  suitable  for  numerical  computation. 

The  propagation  of  a stress  wave  of  known  characteristics 
through  a cylindrical  specimen  held  between  the  flat  ends  of  two 
elastic  bars  is  analyzed.  Three  specimens  are  considered,  namely, 
one  made  of  a purely  elastic  material,  and  a second  and  third 
made  of  materials  which  display  elasto- viscous  and  hereditary 
characteristics,  respectively.  An  approximate  solution  of  the 
problem  is  given  in  the  report  for  each  of  these  specimens. 
Further,  a more  elaborate  and  complete  theory  of  the  same  problems 
is  indicated.  In  Appendix  A a mathematical  analysis  of  the 
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longitudinal  impact  of  a ball  on  an  inf initely^long  cylindrical 
bar  is  given „ The  final  results  of  this  analysis  are  presented 
in  tabular  form  for  given  initial  velocities  of  the  impinging 
balls. 
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AT  LOW  RATES  OF  COMPRESSIVE  STRAINING 


SECTION  II 

Description  and  Use  of  the  Apparatus  for  Testing  Materials 

at  Low  Rates  of  Straining 

A.  General  description  of  the  apparatus 

The  apparatus  employed  in  these  experiments  is  shown  in 
Figures  1 amd  2 and  in  the  schematic  diagram  of  Figures  3 and  4, 
It  consists  essentially  of: 

1)  two  3-foot,  steel  bars  of  equal  mass  suspended  as 
ballistic  pendulums; 

2)  a rotating-drum  camera,  the  drum  of  which  rotates  at 
a known  speed,  and  the  shutter  of  which  is  synchronized  to 
operate  with  the  motion  of  the  steel  bars; 

3)  an  optical  system  which  focuses  the  image  of  a very 
thin  slit  on  the  knife  edges  machined  on  the  adjoining  ends 
of  the  steel  bars  (a  cylindrical  specimen  to  be  tested  is 
mounted  between  the  two  bar  ends  on  which  the  knife  ends  are 
machined) ; 

4)  an  electromagnetic  device  which  can  r elease  one  or 
both  of  the  bars  at  the  same  time. 

The  cylindrical  specimens  of  the  rubber  or  rubber-like 
materials  are  1/2-inch  in  diameter  and  1/2-inch  long.  The 
first  group  of  specimens  has  been  fabricated  from  the  following 
materials: 
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Fig.  2 BACK  VIEW  - CAMERA 


SCHEMATIC  DIAGRAM  OF  EQUIPMENT  ARRANGEMENT  - OPTICAL  METHOD 


Figure  4.  Schematic  diagram  of  optical  system 


P.  . o ..Projection  Lamp,  GE  750T12P 

Go 0 0 .Ground-Glass  Plate,  2 in.  x 3 in. 

S.„..Slit,  0.001  in.  x 7/8  in. 

Ko.o. Knife  Edges  on  bars 

M..„. Mirror,  first-surface  set  at  45° 
to  axis  of  lens  arrangement 

F. a. .Film,  Kodak  Linagraph  Panchromatic  LP421 

L-p  o .Condensing  Lens,  6 in.  f/1,  4-7/16  in„  diam. 

Lgc o 0Projection  Lens,  75mm  f/1,  49mm  diam  coated 
achromat 

LgoooPlano  Convex  Field  Lens,  17  in.  f/1,  3-1/2  in. 
diam. 

L^, •• Plano  Convex  Field  Lens,  21  in  f/1,  3-1/2  in. 
diam. 

L^0 o .Camera  Lens,  Color  Skopar  f/3.5,  105  mm  f/1 


Side  View 


Fig.  4 SCHEMATIC  DIAGRAM  OF  OPTICAL  SYSTEM  - OPTICAL  METHOD 
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a.  Natural  rubber  (27-30  Durometer  hardness) 

b.  Butyl-25  (35-40  Durometer  hardness) 

c.  Neoprene  GN  (35-40  Durometer  hardness) 

d.  Hycar  OR-15  (35-40  Durometer  hardness) 

e.  Buna-N  (35-40  Durometer  hardness) 

f.  GRS-1585  (35-40  Durometer  hardness) 

The  recipes  for  the  rubber  specimens  are  given  in  Appendix  CL. 
The  second  group  will  be  fabricated  from  the  same  materials, 
but  in  about  70  Durometer  hardness. 

The  specimens  are  placed  on  the  plane  end  of  one  of  the 
steel  bars  such  that  the  longitudinal  axis  of  the  bar  and  the 
specimen  coincide.  The  other  steel  bar  is  released  from  a 
predetermined  height  and  made  to  impinge  upon  the  free  end  of 
the  specimen.  During  the  impact,  a photograph  is  taken  of  the 
interval  between  the  knife  edges  which  lie  in  the  plane  of  the 
ends  of  the  steel  bars.  The  film  is  calibrated  by  taking  a 
still  photograph  of  the  interval  when  both  bars  are  barely  in 
contact  with  the  specimen. 

The  deformation  of  the  specimen  at  succeeding  instants  of 
time  is  measured  directly  on  the  film  by  means  of  a traveling 
microscope.  This  instrument  is  fitted  with  two  orthogonal, 
independent  motions,  each  calibrated  to  read  a total  excursion 
of  4,5  inches  in  increments  of  0,001  inch.  From  these  measure- 
ments it  is  possible  to  deduce  the  dynamic  stress-strain 
relationship  of  the  specimen. 

With  certain  modifications,  the  same  apparatus  will  be 
used  to  test  specimens  at  ambient  temperatures  other  than  room 
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temperature o A small  chamber  will  be  constructed  which  will 
completely  enclose  the  specimen  and  yet  not  interfere  with 
any  of  the  moving  parts  necessary  to  conduct  the  tests „ Pre- 
liminary  designs  of  the  chamber  have  been  considered  and  it 
will  be  constructed  when  the  testing  program  at  room  tem- 
perature has  been  completed,,  With  acetone  and  dry  ice  (or 
possibly  another  solvent)  it  should  be  possible  to  conduct 
tests  at  ambient  temperatures  of  -68F,  and,  with  hot  water, 
at  +158F.  These  devices  are  not  described  in  this  report  but 
will  be  discussed  more  completely  in  a later  technical  report „ 

A static  testing  device  has  been  built  and  is  shown  in 
Figure  5.  Static  loads  up  to  200  pounds  can  be  applied  to  the 
specimen  and  a maximum  deformation  of  1-inch  can  be  measured 
to  within  0.00025  inch.  The  same  size  specimens  (1/2-inch  x 
1/ 2-inch  cylinders)  are  used  for  the  static  tests.  Compression 
set  tests  of  each  rubber  stock  are  being  conducted  concurrently 
by  the  staff  of  the  Rubber  Laboratory,  Armour  Research  Founda- 
tion, an  affiliate  of  IIT. 

B.  Theory  of  the  Method 

In  Figure  6 let  F(t)  represent  the  compressive  force  acting 
on  the  specimen  at  the  instant  t measured  from  the  instant 
t = 0,  at  which  time  the  impinging  bar  first  contacts  the  speci- 
men. Let  x-,  (t)  and  x2(t)  be  the  horizontal  displacements  of 
the  two  extremities  of  the  bars  measured  from  their  position  at 
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the  instant  t * 0 and  m be  the  mass  of  each  of  the  two  equal 
steel  bars. 

The  following  relationship  will  be  verified  [5],  [6],  [7], 


F * - m 


d (x2  - 

— 5* 


xl> 


m 


d2s 


(2.1) 


Gall  4(t)  the  actual  length  of  the  specimen  at  the  instant  t 
U is  deduced  directly  from  the  readings  of  the  film) . 

Then 


l (t)  * t>  Q - s (t) 

where  l is  the  initial  length  of  the  specimen. 

By  differentiating  the  experimental  ( s,  t)  curves  twice  it 
is  possible  to  find  the  corresponding  ( F,  t)  curves  by  the  use 
of  equation  (2.1).  If  we  suppose  that  during  deformation  the 
volume  of  the  specimen  remains  constant  and  call  A and  Aq 
the  cross-section  areas  of  the  specimen  at  the  instants  t and 
t * 0,  respectively,  and  VQ  initial  volume,  we  get 


A * = A* 

o o 


The  actual  true  stress  will  then  be 


a(t) 


F (t)  _ F (t ) *(t) 

KTET 


as  compared  to  the  engineering  stress 


»(t)  -m 


o 
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The  actual  true  strain  is  defined  as 
£ (t)  = In 

as  compared  to  the  engineering  strain 

e (t)  - SW 
x o 

From  the  readings  of  the  photographs  it  is  possible  to 
determine  the  curves:  t (t)  „ F(t)a  cr(t),  e(t)  and  finally  the 

dynamic  stress-strain  curves  (a,  e) . 

C.  Degree  of  aqcuracy  of  the  mechanical  and  optical  methods 
employed 

Two  causes  of  possible  error  in  the  described  experimental 
method  have  to  be  taken  into  account*  The  first,  of  a mechanical 
nature , is  due  to  the  elastic  deformation  of  the  steel  bars 
during  impact.  The  second,  of  an  optical  nature,  is  due  to  the 

finite  dimensions  of  the  slit* 

The  effect  of  these  causes  of  error  are  negligible  in  the 
case  of  plastics  and  rubber-like  materials  due  to  the  fact  that 
during  impact  the  forces  acting  on  the  metallic  bars  and  the 
induced  deformations  of  the  bars  are  very  small,  while  the  Lotal 
duration  of  impact  is  comparatively  large.  An  evaluation  of  their 

effects  is  given  below. 

1)  Let  the  following  notation  be  used; 

the  density  of  steel 

E,  Young" s modulus  of  steel 
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velocity  of  sound  in  steel 

pressure  of  the  longitudinal  wave  produced 
in  bars  by  the  force  F(t)  caused  by  the 
impact 

the  particle  velocity  of  the  same  wave 

acceleration  of  the  bar  at  the  instant  of 
initial  impact  with  the  specimen 
velocity  of  the  bar  at  the  instant  of 
initial  impact  with  the  specimen 
cross-section  area  of  the  bars 


(2.2) 


(2.3) 


(2.4) 
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Let  vQ  be  the  velocity  (negative)  of  the  bar  after  impact.  Then 

from  equation  (2.4)  the  maximum  elastic  deformation  in  the  bar 
is  given  by: 

L^  , 2 Lv 

5max  " 5^  Cto  + To  3 (2.5) 

which  yields  an  approximate  value  for  \ of  0.01-inch  for 

IQclX 

the  bars  used. 

2)  The  finite  thickness  of  the  slit  image  of  the  film 

gives  rise  to  an  error  in  determining  the  instants  of  time  for 

the  experimental  s(t)  curves.  The  error  is  of  the  same  order  of 

magnitude  as  the  time  taken  by  the  film  on  the  rotating  drum  to 

move  a distance  equal  to  the  thickness  of  the  slit  image  on  the 

film  and  can  be  determined  as  follows: 

Call  8 the  thickness  of  the  image  of  the  slit  on  the  film, 

L,  the  total  length  of  the  film,  n the  number  of  revolutions  per 

second  of  the  rotating  drum.  It  follows  that  a 1"  length  of 

film  corresponds  to ^ seconds  and  the  thickness  of  the  slit 

* 

image  corresponds  to  ^ seconds.  If  T is  the  total  duration  of 

impact,  the  percentage  error  in  the  evaluation  of  time  in  our 

ft 

measurements  will  be  given  by  which  indicates  a possible 

error  in  the  indicated  time  of  approximately  0,2  percent. 

D.  Direct  computation  of  the  dynamic  stress-strain  curves  for  a 
given  material,  with  known  hereditary  characteristics,  under 
impact  loading. 

When  the  hereditary  characteristics  for  a given  material  are 
known,  as  a function  of  time,  it  is  possible  to  calculate  directly 
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the  dynamic  stress-strain  curve  under  impact . In  order  to 
simplify  the  calculations  we  shall  assume  that  under  static 
load  conditions  the  material  obeys  Hookers  law,  while  under 
dynamic  conditions  the  stress-strain  relationship  for  the 
material  is 

c = Ee  + *(t  - t)  dr  (2.6) 

Suppose  that  such  a material  is  subjected  to  an  impact 
test.  Call  £ its  initial  lengthy  Aq  its  initial  cross-sectional 
area,  £(t)  ~l  Q “ s(t),  its  length  at  time  t 0 and  m the  mass 
of  the  impinging  bar. 

The  following  integro-differential  equation  will  characterize 
the  motion  during  impact  in  the  case  of  the  assumed  linear  stress- 
strain  relationships 

nt  - x)  dT  = 0 

dt2  lo  lo  0 

(2.7) 

The  solutions  of  equation  (2.7)  will  be  given  for  different  types 
of  heredity.  It  is  evident  that  if  the  equation  is  solved,  i.e., 
the  function  s(t)  is  determined,  then  the  stress-strain  curves 
can  be  calculated  by  the  use  of  the  formulas  given  in  Section  II  B. 

1.  Heredity  of  the  first  degrees  V (t)  = He 

By  a simple  differentiation  with  respect  to  t equation  (2.7) 
can  be  transformed  into  the  following  differential  equation; 
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m 


d3s (t ) . _ _ d2s(t)  ^ (E  + H)  Aods.  AoBas(t)  _ A 


dtv 


dtJ 


to 


£ o 


(2.8) 


The  solution  of  equation  (2.8)  is  of  the  form 
s(t)  = Ci  e”nt  + C.,e"ptcosqt  + C,>e~ptsinqt 


(2.9) 


where  -n;  "P + iqj  -p-iq,  are  the  roots  of  the  algebraic  equations 


r3  + ax2  + x + 5#lJS.  . o 


7 m 

^ o 


In* 


(2.10) 


From  the  initial  conditions? 


s = 0,  gf  = * 0 for  t * 0.  it  follows  that 

dt 


Gi  =Ylvo‘  C2  = “yTvo*  C3  =y2v 


lvo*  3 r 2 o 


(2.11) 


where  s 


2p 

2 77  IT 

q + (p-n) 


-2  * n2 


q(qZ  + (p=n)  ) 


(2.12) 


2.  Heredity  of  the  first  degree  with  residuals 
t (t)  = He~at  + H0 

The  solution  of  equation  (2.8)  is  given  by  equations  (2.10) # 
(2.11) , and  (2.12)  where,  in  this  case,  =n,  -p* iq,  “P  -iq  ar«  the 
roots  of  the  algebraic  equations 
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« 9 (E  + H ♦ H )A  A (E  ♦ H )a 

O . Z , 0 0*0  o 

X + x + — — -r-— X + — — i-'TT 

£ om  £ om 


(2.13) 


3.  Heredity  of  the  second  degree:  ^(t)  = H-,  e al^  + H9e  a2c 

By  a double  differentiation  with  respect  to  t equation 
(2.8)  can  be  transformed,,  in  this  case,  into  the  following 
differential  equation: 


H + w x . „ „ -id2s (t) 
Hi  * ^2  1 dt^  " 

EA 

+ rn  tE(ai + a2>  + “A  + °2HP  It + aia2  Hr  ■ 0 (2-14) 


dMti  * ( ♦ , Aw  + ^ 

^ 1 A dt6  rom 


dt 

Ao 

mxzz 

O 


The  solution  of  equation  (2.14)  is: 

s ( t ) - G je~  111  + C2e'"’^t  + C3e”ptcosqt  + C4e”ptsinqt 

(2.15) 


where  - t),  - C , =p  + iq„  -p  - iq  are  the  roots  of  the  algebraic 
equation: 


>)x' 


(E  + Hx  + H2) 


+ 


ala2 


+ a2)  + a1H2 


+ 


E A. 


a2Hx]x  + aia2  r 


— - 0 


m 


(2.16) 


The  initial  conditions  for  determining  the  four  constants  C^,  C^c 
G^  and  C4  ares 
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S - Os  - vQ°a 


d3s  <E  + H1  + H2)Ao  ds 


for  t * 0 

The  following  four  equations  for  determining  the  constants  Gj 
(i  = 1„  2,  30  4)  are  derived? 


C1  * C2  * C3 


■riC-,  - pCo  + qCA  = v 


rl2C1  + C 2C2  + (p2“q2)G3  - 2pqG4  - 0 


( r)°  +r)B)C1  - I C°  + CB)C2  + p(3q2°p2°B)C3  + q (3pi5-qa+B)C4  - 0 


(2.17) 


2 2, 


where? 

(E  + + H2)Ao 

B - - — ~ 

4c  Heredity  of  the  second  degree  with  residual r 

*(t)  - H-^^l1  + H2e“a2t  + H0 

The  solution  of  equation  ( 2 D 8 ) is  given  by  equation  (2„15) 
and  (20 17)  where?  -t),  - £ „ ~p  + igff  -p  ~iq  are  the  roots  of  the 
algebraic  equations 


!“1  + a2)x3  + !E  + H1  + H2  ♦ H0)  -♦  aia2]x2 


(2.18) 


t (E  - H0>  (ax  ♦ a2)  +.0^  - o^]*  + (E IJslh.  = 0 
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and  B of  equation  ( 2 « 1 7 ) is  given  by: 

R . (E  + Hi  + H2  + Ho)A0 
— 

E . Determination  of  the  hereditary  characteristics  of  a material 
from  the  experimental  dynamic  stress-strain  curves 

We  shall  demonstrate  the  method  of  analysis  by  considering 
the  case  of  heredity  of  the  first  degree  with  residual. 

By  equation  (2.9)  we  note  that  s(t)  can  be  written  in 
the  form 

s ( t ) = G1e”nt  + Ce“Pt  + U e (2.19) 

where  and  n are  real,  C and  P are  complex,  and 

G - (~C1  * iC3}  , p - p ♦ iq  (2.20) 

2 

We  assume  an  expression  for  s(t)  of  the  form  of  Equation  (2.19) 
and  fix.  it  to  our  experimental  data  by  the  method  of  least  squares, 
following  the  procedure  of  Prony  (See  reference  [8],  page  373). 
Since  -n,  -p„  and  -F  are  the  roots  of  the  equation 

x3  + (n  + 2p)x2  + (2pn  + p2  ♦ q2)x  + n(p2  + q2)  * 0 (2.21) 

We  have  by  comparison  of  (2.13)  and  (2.21)  the  following  three 
equations  for  determination  of  a t HQ,  and  H: 

a = n + 2p 

2 2 

E + Hq  + H - (2pn  + p * q ) 
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“•tn  9 9 

a(Hc  + E)  - -fS  (pz  + qz)n 

where  E is  the  static  secant  modulus  of  the  material. 

Since 

af^t-o  = vo  = qC3  ‘ (n“p)Gi 

we  have  a check  on  the  accuracy  of  our  calculations,  as  we 
can  obtain  the  initial  velocity  directly  from  the  experimental 
data. 

With  obvious  modif icati-ons  this  method  can  also  be  used 
to  analyze  viscous  damping. 

However,  it  should  be  noted  that,  in  the  relationship 

, = De  + 

> 

D is  a lumped  parameter  which  is  obtained  directly  from  the 
dynamic  test  data?  it  may  have  no  relationship  to  the  static 
secant  modulus  E„  A close  approximation  of  its  value  is 
D = E + Hq. 
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DESCRIPTION  AND  USE  OF  APPARATUS  FOR  TESTING  MATERIALS 


AT  HIGH  RATES  OF  STRAINING 


SECTION  III 

Description  and  Use  of  Apparatus  for  Testing  Materials 
At  High  Rates  of  Compressive  Straining 

A,  Ge Mral  description  of  the  mechanical  apparatus. 

The  mechanical  apparatus  is  shown  in  Figure  7.  It 
consists  of: 

1)  two  steel  bars,  1"  diameter,  6 feet  long,  and 

2) .  a short  impinging  bar  (hammer)  with  a round  head. 

This  mechanical  apparatus  is  equipped  with  a system  of  strain 
gages,  amplifier,  radio  frequency  oscillator,  oscilloscope  and 
photographic  camera  to  record  the  longitudinal  pulses  which  are 
propagated  along  the  two  steel  bars  when  the  short  bar  strikes 
the  end  of  one  of  the  long  bars.  The  details  of  the  electronic 
equipment  are  given  in  the  next  section. 

The  specimens  to  be  tested,  in  the  form  of  circular  discs 
of  uniform  thickness,  with  diamefers  slightly  smaller  than  those 
of  the  adjoining  bars  are  placed  between  the  carefully-lapped 
plane  ends  of  the  two  bars.  For  reducing  friction  between  speci- 
men and  bars  both  surfaces  of  each  specimen  are  lubricated  be- 
fore being  placed  in  the  apparatus. 

The  hammer  is  released  from  a known  height.  When  it  strikes 
the  first  steel  bar  an  open  ground  circuit  is  completed  which 
triggers  one  sweep  of  the  dual-beam  oscilloscope.  A compression 
wave  is  propagated  in  the  bar  and  its  passage  is  recorded  from 
an  amplified  strain  gage  signal  by  the  camera  attached  to  the 


26 


VIEW  OF  ASSEMBLED  EQUIPMENT  - STRAIN-GAGE  METHOD 


- 28 


oscilloscope „ When  the  compression  wave  reaches  the  other  end 
of  the  steel  bar,  a certain  portion  of  the  wave  is  reflected 
back  along  the  -bar  as  a tension  wave  and  the  remainder  is 
transmitted  through  the  specimen  on  to  the  second  steel  bar„ 

On  the  first  steel  bar,  about  12  inches  from  the  mounted 
specimen,,  a strain  gage  is  fastened,,  It  is  connected  electri- 
cally  to  a second  trigger  circuit  which  sets  the  second  sweep 
of  the  dual-beam  osci lloscope 0 Thus,  when  the  compressive  wave, 
which  is  being  propagated  in  the  second  steel  bar,  reaches  a 
set  of  strain  gages  about  12  inches  beyond  the  specimen  a second 
compressive  trace  is  recorded  by  the  oscilloscope  cameraD 

From  this  photographic  record  then  it  is  possible  to  cal- 
culate  directly  the  dynamic  stress-strain  curve  of  the  specimen, 
as  will  be  shown  later  „ 

In  this  method,  auxiliary  electronic  equipment  is  used  to 
amplify  a signal  from  a set  of  strain  gages,  to  superpose  a 
timing  pulse  on  this  signal,  and  finally,  to  indicate  the  mag- 
nitude and  length  of  the  impact  pulses  on  the  screen  of  an 
oscilloscope „ 

Essentially,  there  are  three  circuits  of  importance, 

1)  the  strain  gage  circuit,  which  amplifies  the  original 
signal  and  mixes  it  with  the  timing  pulse  before  it  is  fed  into 
the  signal  circuit  of  one  beam  of  the  dual-beam  oscilloscope? 

2)  the  time  pulse  circuit,  which  superposes  20  microsecond 
pips  on  the  strain  gage  signal?  and 
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3)  the  trigger  circuits  which  set  the  sweep  of  one  beam  of 
the  oscilloscope  to  fire  when  they  receive  an  external  voltage 
pulse. 

The  strain  gages  are  arranged  in  sets  of  four  at  equi- 
distant points  around  the  periphery  of  the  steel  bars  (Figure  7 ) . 
There  are  two  sets,  each  about  12  inches  from  the  head  end  of 
the  two  steel  barsD  Since  each  set  is  connected  in  series,,  the 
net  effect  of  any  bending  modes  will  be  zero  and  the  voltage 
change  generated  will  be  due  entirely  to  the  compressional  pulse 
being  propagated  along  the  steel  bars. 

The  voltage  pulse  is  amplified  first  by  a four-stage  triode 
amplifier  with  a wide  frequency  response „ A timing  pip  at  in- 
tervals of  20  microseconds  is  added  to  the  output  of  the  ampli- 
fier. 

The  time  interval  of  the  pip  is  controlled  by  a Hewlett- 
Packard  audio-oscillator  and  can  be  varied  if  it  is  desired. 

The  pulse  (pip)  generator  consists  of  an  884  gas  triode  and  a 
differentiating  circuit.  The  output  of  the  latter  is  formed  by 
a clipper  circuit  and  then  fed  into  a wide  range  amplifier  in 
order  to  obtain  the  proper  phase  relationship  between  the  original 
pulse  and  the  timing  pip.  This  corrected  pulse  is  superposed, 
in  turn,  on  the  amplified  strain  gage  signal  just  before  the 
latter  is  fed  into  one  of  the  signal  circuits  of  the  oscillo- 
scope. 

The  phenomenon  is  transient  rather  than  steady  state  and  as 
a consequence  it  is  necessary  to  trigger  the  sweep  of  the  oscillo- 
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scope  just  before  the  aompressional  pulse  reaches  the  strain 
gages.  This  is  accomplished  by  making  use  of  the  circuits  within 
the  oscilloscope  and  additional  external  circuits. 

The  sweep  oscillator  for  each  beam  of  the  oscilloscope  is 
set  for  a pre-determined  sweep  time,  and  then  the  other  controls 
are  so  adjusted  that  neither  oscillator  will  fire  until  an 
additional  external  voltage  pulse  is  impressed  upon  them.  Two 
external  trigger  circuits  are  used  to  insure  that  the  oscillators 
will  fire  at  the  proper  time. 

The  circuit  for  the  first  bar  is  energized  when  the  hammer 
completes  a ground  circuit  at  the  instant  of  initial  impact;  for 
the  second,  when  the  compressional  wave  passes  a strain  gage  about 
20  inches  in  front  of  the  second  set  of  recording  strain  gages* 

In  either  case,  a voltage  pulse  is  amplified  and  then  applied  to 
a sweep  oscillator  of  the  oscilloscope  in  order  to  produce  a 
single  sweep  at  the  instant  the  compressional  wave  passes  the 
recording  strain  gages.  Thus,  only  one  pulse  from  each  set  of 

i 

strain  gages  is  recorded  by  the  camera. 

This  method  of  testing  will  allow  determination  of  the 
stress-strain  curves  for  rates  of  straining  of  the  order  of  micro- 
seconds. The  length  of  the  direct  pulse  in  the  first  long  bar 
is  governed  by  the  length  of  the  striking  bar  (hammer)  and  the 
relative  magnitude  of  the  pulse  by  the  velocity  at  the  time  of 
initial  impact* 
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In  the  following  table  we  give  some  of  the  computed  character 
istics  of  the  pulses,, 


TABLE  I 

Theoretical  time  of  contact  in  microseconds  as 
a function  of  height  of  drop  and  hammer  length 


Drop  height 
Hammer  length 

5 cm 

10  cm 

15  cm 

25  cm 

25  cm 

140 

136 

134 

130 

50  cm 

239 

233 

231 

228 

75  cm 

336 

331 

328 

325 

100  cm 

432 

427 

424 

421 

experimental  method  of 

determining  dynamic  stress-strain 

curves  for  very  high  rates  of  straining  is  subject  to  certain 
restrictions,  as  is  the  method  based  on  the  use  of  the  Davies- 
Hopkinson  pressure  bar  [10],  [11],  namely „ 

1)  The  waves  propagated  in  the  bar  are  assumed  to  be  elastic 
waves,  i „e „ , the  stress  at  every  point  in  the  bar  must  always  be 
within  the  region  where  the  stress-strain  curve  is  linear  and 
reversible,, 

2)  It  is  assumed  that  the  pressure  wave  is  propagated 
without  distortion „ This  assumption  is  only  true  when  the  wave 
lengths  of  the  elastic  waves  are  large  compared  with  the  lateral 
dimensions  of  the  bars#  When  this  condition  is  not  fulfilled. 
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the  wave  suffers  dispersion  and  the  form  of  the  pulse  is  dis- 
torted as  it  travels  along  the  bar„ 

3)  A third  assumption  implicit  in  the  method  is  that  the 
pressure  in  the  pulse  is  uniformly  distributed  over  the  cross- 
section  of  the-bar,  even  when  the  force  acting  on  the  end  is 
concentrated  over  a small  area  surrounding  the  center,, 

However,  this  method  of  investigation  presents  the  follow3 
ing  advantages  in  comparison  with  the  method  using  the  Davies- 
Hopkinson  Bars 

1)  While  in  the  Davies  method  it  is  possible  to  record 
only  the  first  part  of  the  stress-strain  curve,  i„e„,  the  part 
when  the  specimen  is  compressed  but  not  when  it  recoils,  with 
this  method  the  complete  stress-strain  curve  can  be  investigated. 

2)  The  determination  of  the  stress- strain  curves  from  the 
experimental  data  is  much  simpler,  and  there  is  less  possibility 
of  making  mistakes  in  the  numerical  computation.  In  fact, 
while  in  the  Davies'7  method  a numerical  differentiation  has  to 
be  performed,  in  this  method  the  stress-strain  curves  are  ob- 
tained through  numerical  integration  of  the  experimental  data, 

a method  of  computation  with  many  inherent  advantages. 

B„  Theory  of  the  method 

When  the  hammer  hits  the  first  bar  a direct  compressive 
wave  aD  (Figure  8)  is  transmitted  undistorted  along  the  first 
long  steel  bar  with  velocity  c2  - = [12].  When  it  reaches 

the  specimen,  after  a time  t = ~ (L  being  the  length  of  the  bar) 
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the  wave  suffers  dispersion  and  the  form  of  the  pulse  is  dis- 
torted as  it  travels  along  the  bar0 

3)  A third  assumption  implicit  in  the  method  is  that  the 
pressure  in  the  pulse  is  uniformly  distributed  over  the  cross- 
section  of  the-bar*  even  when  the  force  acting  on  the  end  is 
concentrated  over  a small  area  surrounding  the  center,, 

However,,  this  method  of  investigation  presents  the  follow* 
ing  advantages  in  comparison  with  the  method  using  the  Davies- 
Hopkinson  Bars 

1)  While  in  the  Davies  method  it  is  possible  to  record 
only  the  first  part  of  the  stress-strain  curve * iDe„*  the  part 
when  the  specimen  is  compressed  but  not  when  it  recoils,  with 
this  method  the  complete  stress-strain  curve  can  be  investigated. 

2)  The  determination  of  the  stress-strain  curves  from  the 
experimental  data  is  much  simpler * and  there  is  less  possibility 
of  making  mistakes  in  the  numerical  computation.  In  fact* 
while  in  the  Davies'7  method  a numerical  differentiation  has  to 
be  performed*  in  this  method  the  stress-strain  curves  are  ob- 
tained through  numerical  integration  of  the  experimental  data* 

a method  of  computation  with  many  inherent  advantages. 

B»  Theory  of  the  method 

When  the  hammer  hits  the  first  bar  a direct  compressive 
wave  aD  (Figure  8)  is  transmitted  undistorted  along  the  first 
long  steel  bar  with  velocity  c2  - | [12].  When  it  reaches 

r 

the  specimen*  after  a time  t = ~ (L  being  the  length  of  the  bar) 


the  wave  suffers  dispersion  and  the  form  of  the  pulse  is  dis- 
torted as  it  travels  along  the  barD 

3)  A third  assumption  implicit  in  the  method  is  that  the 
pressure  in  the  pulse  is  uniformly  distributed  over  the  cross- 
section  of  the- bar,,  even  when  the  force  acting  on  the  end  is 
concentrated  over  a small  area  surrounding  the  center,, 

However,  this  method  of  investigation  presents  the  follow3* 
ing  advantages  in  comparison  with  the  method  using  the  Davies- 
Hopkinson  Bars 
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of  ma^Aiiy  miataices  in  tne  numerical  computation.  In  fact, 
while  in  the  Davies'7  method  a numerical  differentiation  has  to 
be  performed,  in  this  method  the  stress-strain  curves  are  ob- 
tained through  numerical  integration  of  the  experimental  data, 
a method  of  computation  with  many  inherent  advantages. 


B . Theory  of  the  method 

When  the  hammer  hits  the  first  bar  a direct  compressive 
wave  (Figure  8)  is  transmitted  undistorted  along  the  first 
long  steel  bar  with  velocity  c"'  38  = [12],  When  it  reaches 

the  specimen,  after  a time  t = ==  (L  being  the  length  of  the  bar) 
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Pulse-Time  Relations  for  Figure 
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only  <Xp,  a part  of  the  original  pulse,  is  transmitted  through  the 
specimen  and  into  the  second  bar.  A tension  pulse  is  reflected 
back  along  the  first  bar- 

The  theory  of  impact  developed  here  can  be  made  more  realistic 
by  considering  a cylindrical  specimen  of  a material  held  between 
the  plane  ends  of  two  steel  bars  (Figure  8) . A stress  wave  of 
known  properties  is  propagated  along  the  first  steel  bar,  and 
the  problem  is  to  determine  the  characteristic  stress  waves  in- 
duced  in  the  specimens  which  follow  various  stress-strain  relation- 
ships. 

A minimum  number  of  simplifying  assumptions  will  be  made  for 
both  a relatively  realistic  analysis,  and  later,  for  a more  appro- 
ximate theory. 

These  initial  assumptions  for  the  first  analysis  are; 

(1)  The  stress  in  the  specimen  is  distributed  uniformly 
across  each  cross  section. 

(2)  There  are  no  reflections  in  the  specimen. 

(3)  aR  (t ) = aD(t)  - qo (o,  t)  for  x * 0 and  all  values  of  t? 

aT(t)  = qo(£,  t)  for  x = l and  all  values  of  t,  where  the 
specimen  is  of  length  £ and  cross-sectional  area  A2,  the  steel 
bar  has  area  A^  and  q — A2/A1  is  the  ratio  of  the  areas.  a^(t) 
is  the  wave  transmitted  initially  along  the  bar  on  the  left;; 
aD(t)  is  the  wave  reflected  from  the  interface  of  the  left  bar 
and  the  specimen?  and  aT(t)  is  the  wave  transmitted  to  the  steel 


bar  on  the  right. 


(2a) 

(Hooke) 
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Figure  9 

Schematic  Diagram  of  Material  Models 


The  third  assumption  is,  essentially,  a consequence  of  the 
first  two  assumptions  and  the  equilibrium  of  the  specimen „ The 
second  of  our  assumptions  allows  us  to  consider  the  specimen 
infinitely  extended  to  the  right 0 

The  stress-strain  relationships  for  the  specimen  materials 
considered  are 

(1)  Elastic  (Hookean  Model) 

(2)  Elasto-Viscous  (Voigt  Model)  [14] 

(3)  Hereditary  (Boltzmann  Model)  [15],  [16] 

and  the  models  which  are  characterized  by  each  of  these  relation- 
ships  are  shown  in  Figure  9.  In  all  cases  the  equation  of  motion 
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is 


-~^x  <i(x,  t)  53  P g u(x,  t) 


where  pg  is  the  mass  density  of  the  specimen, 


(3.1) 


Boundary  and  Initial  Conditions 

The  basic  boundary  and  wave  front  conditions  are  the  same 
for  all  three  of  the  stress-strain  relations  considered,  namely. 


For  t ^ 0;  0 ^ x^  oo  a “ 0 


(3.2) 


For  t ^ 0 the  particle  velocities  of  the  specimen  and  the  steel 
rod  are  the  same,  that  is  for  t ^ 0,  x = 0 

W = V1  " C°D(t)  + aR{t)l  (3,3) 


where  u(x,  t)  is  the  displacement  of  the  specimen;  and  c^ 
are  the  density  and  wave  velocity  of  the  steel  rods,  respectively. 
The  known  initial  wave  is  assumed  to  have  the  form 


aD (t)  = 


0 for  t ^ 0 
b sin\t  for  0 ^ t ^ T 


L 0 for  T ^ t 


(3.4) 


where 

■x  _ nn 

x ' T 

n is  an  integer  and  T is  the  time  length  of  the  pulse.  The 
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linearity  of  our  equations  permits  superposition  of  solutions  for 
different  n. 

Case  I - Elastic  Specimen 


a(x,  t)  53  E2  -gj  u(x„  t) 


(3.5) 


where  Eg  is  the  static  Young" s modulus  for  the  specimen  material, 
Elimination  of  u(x,  t)  between  (3.1),  (3.5)  yields 


1? 


2 d2a 


(3.6) 


The  first  assumption  together  with  conditions  (3.2),  (3,3)  and 
equation  (3.5)  yield  the  boundary  conditions. 

Our  system  becomes 

a2a  _ 2 a2a 

atz  z ax^ 


t - -5-*o 
c2 

0 ^ x 


a = 0 


t3.7) 


t > 0 


x = 0 


a(o,  t)  = 


| °D(t) 


9 p2c2 
P 2°  2 

^lcl 
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In  view  of  our  assumption  of  no  reflections  in  the  specimen,  the 
solution  of  (3,7)  is 


a(x,  t)  * 0 for  t - 


x < 


0,  0 $ x 


a(x,  t)  * ~ an(t  - ~)  for  t ^ 0,  0 < x <£ 
q u c ^ 


(3.8) 


Since  the  transmitted  stress  a rp  is  given  by  cr,j,(t)  = qcr(£,  t) 
we  have  the  transmitted  wave  as  an  immediate  consequence  of  (3.8). 
For  the  initial  wave  given  in  (3,4)  we  have 

a>p(t)  = 0 for  t “ - 0 1 


Om(t)  = 0bsin\(t  - ) for  0 ^ t ^ T 

i c2 


(3.9) 


aT(t)  * 0 


for  T - t 


Case  II:  Elasto-Viscous  Specimen 

du(x,  t)  <32u(x,  t) 

<j(x,  t)  - D2  <3x  h Y 3t'3x 


(3.10) 


where  Y is  the  coefficient  of  viscosity.  The  initial  conditions 


are : 


t =*  0:  u(x„  o)  - °)  “ = ^ 

, v da(x„  o)  _ n 
<y(x,  o)  - ~ 0 


(3.11) 


The  boundary  conditions  are: 


= 0:  Jr  u(o,  t)  = “p"V*  C 2od (t ) -qcr(o,  t)]  (3.12) 
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where  an(t)  is  given  by  (3.4) 


X— prCOi  u(x,  t ) — >-0 

aUs  t)->0 


(3.12) 

cont. 


We  attack  this  problem  by  means  of  the  Laplace  Transform  [17], 
[18],  We  denote  the  transform  of  a sufficiently  regular  function 
f (x,  t)  by  7(x„  z),  i.e.. 


7(x,  z)  « /°°  f(x,  t)  e zt  dt 


f(x,  t)-^7(x,  z) 


(3.13) 


By  Section  4,  [17]  (3.1)  „ (3„10),  (3.11)  „ (3.12),,  imply 
®(x#  z)  = (i>2  y z)^  u(x„  ^ ( 


a(x»  z)  * zzp  2u(x„  z) 


= Os  z u(o„  z)  = - 


pp^2oD  - CI<’(o'  z)] 


(3.14) 


>00  : u(x.  z)—^  0 


a(x„  zh ►O 


where,  by  (3.4) 


aD(z) 


- X>bTl  ° (-l)n  e~zt' 


(3.15) 
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The  solution  of  (3.14)  is 


2p9za,e“<i)ox 
o(x,  z)  - 21 


Vlcl  +p  2zq 


c\, . X 

aSi  a-,  e o 

u(x  z)  = - 0 4 

cooplcl  + p2^ 


(3.16) 


[p2z 


r 


where  coq  is  the  branch  ofit2Wi  which  has  positive  real  part. 


If 


F(x,  t)— : 


D 2-ze  o 


• rc2C^pi+<o0]  <\2  * A 


(3.17) 


then,  by  section  10,  £17] 


F(x,  t)  for  0 = t - T 


t>  s p 2cV  " \ 
P1C1C2 


(3.18) 


F(x,  t)  - (~l)nF (x,  t-T)  for  t - T 


where 


c2  = 


(3.19) 


by  Chapter  VI,  [17] 


D, 


c+iCO 

s 


■’(u  X zt 
ze  ° e dz 


2niVc2  r^P2z 


c-i a>[^4“+a)  J(X2  + *2> 


■p-,c 


F(xf  t) 


(3.20) 
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Let 


g 


D2 

-7  t;  C - 

TT^“  » h “ 


qp2c2 

Plcl  ' 


P 2C2X 


2 C 

(0=1^- 


(3.21) 


then 


F(X.  t)  = 1 / Cep(CT  -^)dC 

2lli  Ti  (<?  + CZ)  (he  + co) 


(3.22; 


where  T ^ is  the  Bromwich  contour  (see  Figure  10)  and  the  branch 
C2 

of  ^ * co  which  has  positive  (or  zero)  real  part  is  chosen. 
This  implies  that  the  only  real  root  of  hC  + to  » 0 for  h > 0 
is  C = 0, 


Case  III:  Hereditary  Specimen 


o (x,  t)  = E 


du  H/1  e"a(t'T)  6 dr 


2 


(3.23) 


*A  method  similar  to  the  above  was  followed  by  I.  N.  Zverev  [19]. 
His  solution  is  given  in  the  form  of  a double  integral  which  must, 
in  general,  be  evaluated  numerically.  For  the  purpose  of  this  re- 
port, numerical  evaluation  of  the  single  integral  (3.22)  is  much 


more  convenient. 
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where  E2  + H is  the  dynamic  Young "s  modulus  and  aB  A are 
Boltzmann  constants  for  the  material.  The  initial  conditions 
are: 


t 0:  u(x,  o)  - (Xj  o)  = ^ (x„  o)  = 0 

0 (x,  o)  = (x0  o)  = 0 


> (3.24) 


The  boundary  conditions  are: 


x = 0:  jp  u(o,  t)  = 


aD(t)  + o(o,  t) 


where  ff^(t)  is  given  by  (3.4) 


oo  s u(xfl  t ) — 

o (xff  t ) ^'-O 

We  take  the  Laplace  transform  of  (30l)p  and  (3.23) 
[17]  and  by  (3.24) , (3.25)  we  have: 


X3.25) 


o(x,  z)  - (E2  + ■£%) 


Hz  ^ au(x^t) 
x 


By  sections  4, 


(a) 


5 — o 

o(x,  z)  = P2z  u(xt  z) 


(b) 


x = 0:  zu(o.  z)  = 


[ 2oD (z)  - qa(o,  z)] 

(c) 


x-^-oo  : u(x„  z)— ►*0 


j(x17  z)— 


(d) 

(e) 


(3.26) 
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and  <Jq(z)  is  given  by  (3.15). 

The  solution  for  (3.1)  and  (3.24)  is 


oU0  z)  * 


2p2poDe^of 


2wae‘“°'i 
u(x,  z)  — «JL_ 


z (P1o1c>o  + qp2^z) 


(3.27) 


where 


2 „ z z ( z + a ) 


E2  + H R _ E2  a 

P2  “ ' Kp+~H  (3.28) 


and  the  branch  of  a)  which  has  positive  real  part  is  taken. 

o 


F(x„  t)-> 


-to. 


+ z / (o)q  + zd) 


(3.29) 


Then,  by  section  10  „ [17] 


a(xff  t)  = 


2b>p2P 

P1G1 


‘ F(x,  t)  for  0 - t — T 


(3.30) 


F(x,  t)  - (“l)nF  (x p t»T)  for  t ^ T 


x *2P  2^ 

where  l * j and  d = j-g- 


By  the  inversion  theorem 
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F (x,  t)  = 


1 


JL 


exp[  zt  - £00  ] z dz 

(7^  + z^)  (a)Q  + zd) 


(3.3lf 


where  ~P 2 the  Bromwich  contour.  We  note,  that,  for  d>0  the 

only  root  of  zd  + coQ  = 0 is  z = 0o  (See  Figure  11) . 

Approximate  Theory 

In  the  following  subsections  we  make  various  physically- 
plausible  assumptions  which  lead  to  solutions  more  amenable  to 
numerical  computation. 

Approximate  Assumptions 

1.  The  stress  in  the  specimen  is  uniformly  distributed 
across  each  cross-section. 

2.  There  are  no  reflections  in  the  specimen. 

3.  The  strain  in  the  specimen  is  uniform. 


4.  The  stress  in  the  specimen  is  independent  of  x,  the 
longitudinal  position  in  the  specimen. 

In  a recent  paper  £ 20] , H»  Charles  presented  a solution  to  the 
problem  of  one-dimensional  hereditary  wave  transmission  which  in- 
volves the  evaluation  of  a fourfold  iterated  integral.  For  the 
purpose  of  this  report,  numerical  evaluation  of  the  single  in- 
tegral (3.31)  is  much  more  convenient. 
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5. 


0R(t)  * (^(t)  - q o ( t ) 


D 


cfm(t)  = q a (t ) 


Case  IV:  Elastic-Specimen  - Approximate 


a (t)  = E2  e 


E2  , 

T (u2 


ti^) 


(3.32) 


by  assumption  3„ 

The  particle  velocity  is  given  by 

V =*  — O 

pc 


(3.33) 


(see  equation  6„  £12]).  Now  the  particle  velocity  of  the  left 
bar  is 

(3.34) 


vi  = ° 9^1  (JC  + aR} 

while  the  particle  velocity  of  the  right  bar  is 


v2 


plcl 


(3.35) 


Thus 


V1  ' HF  (u2 


“l>  = ppj  (oD  + aR  ' °T)  (3>36) 


By  assumption  (5)  we  have 


T2  ' T1  = (°D  ' q0) 


(3.37) 
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Differentiating  (3=32)  and  using  (3=37)  we  obtain 
d 2E2q  2E« 

+ £ ft  ss  ^ 

^ iPlcl  ^Plcl  ^ 

Let 

2E2q 

u 

then  t > 0 implies 

0)  “0)t  cot  i \ , 

o'  = — e S e aD(x)dx 

Now 

J ' e^  sinXxdx  — — -w™” — w£coe^  smX.t"Ke^  co£^»t+\,j 

o co  +X.Z 


Thus 


r 

0 for  t ^ 0 


a (t)  ~ { Y”[cosin\t  - \cos\t  +?v.e”wt]  for  0 ^ 


q (co  + X ) 


[1  “(°l)n  e^1]  e"wt  for  t ^ T 

q (co  + \ ) 


where 


2E2  q 


(3=38) 


(3=39) 


(3=40) 


(3=41) 


t S T 
(3.42) 


to 


(3.43) 
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The  transmitted  stress  is,  thus 
Oj,  “ qcr 


(3.44) 


Case  V:  Elasto-Viscous  Specimen  - Approximate 

'■  D2 

0(t)  - — (u2  - ux)  * * (r2  - vx)  (3.45) 

by  assumption  3.  Differentiate  (3.45)  and  use  (3.37)  obtaining 

da  + 2P2qa  . 2D2°D  + 2 y daD  ,«  46) 

cff  £p  1c1  + 2yq  ^P]_Gi  + 2y  q £P;lc]_  + 2yq  dt 


Let 

2D2q  . ..  _ 2yq 

w i?ici  + 2y^  ' P £p  + 2yq 


(3.47) 


Thus 


qa  (t ) = p,  a^  + (1  - n e f e a^  (t)  dx 


(3.48) 


Hence,  using  (3.41) 


C 0 for  t ^ 0 

— r^sin^-t  + - t1I'4((1)sinU  -^cosM  + X e-^)  ] 
q L a)2  +\2 

a (t)  - ^ for  0 ^ t ^ T 


(3.49) 


b (l^JgK  C1  -(-l)n  •"*]  • "mt 

q co2  + \2 


for  t £ T 
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Case  VI: 


a ( t ) 


Hereditary  Specimen 


- Approximate 
/ e -a(t'T,(v2 


v^)  dr 


(3.50) 


by  assumption  3, 

Eliminate  the  integral  between  (3.50)  and  its  derivative  with 
respect  to  time.  Then  use  (3.37)  obtaining 


d2o  . r . 2q*E2  + H'-,do  ^ 2qaE2  _ 2*E2  + H*  d0D 

* » + zv  xo - — igr  ' Ip  -mr 


+ 


aD 


(3.51) 


Let 


C 


2q(E 

T 


2 * H) 
1C1 


\i  = y (a  + C) 


2E2Hq 

£p?T 


1 (3.52) 

. 


then 


0 for  t ^ 0 


dig,  +2n  da  + (n2  - T]  2)a  - J ~[  (\i2  - h 2)sinXt  + \t  costa] 


— 7 HT 

dtz  ax 


for  0 4 t T 


(3.53) 


v. 


0 


for  t ^ T 
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a ~ Hi’  = 0 for  t ® 0 
o,  continuous  for  t ■ T 


(3.54) 


Thus 

0 for  t ^ 0 

l-b[g  coshrjt  + =■  (^tg  + \h.)  sinhTytle^* 
qo  » oT  * >< 

+b[gcos\t  + h sin'Kt]  for  0 ^ T (3.55) 

be~^  1 -^2™, -JL^L  [ (-l)nelltsinhT)(t“T)  -sinhrjt] 

L T i 

+g[  (' l)ne^  cosh(t“T)  -coshr)t]|-  for  t ^ T 

- w;(ii2  - n2  - X2)  -2uMu2  - n2) 

(\iZ  - r\  - X2)2  -*•  (2\iK 

(3.56) 

h = k -ii  ~>V~-  -’-V  ^ 

((1  - T)2  - X2)2  * (2>u)2 


C,  Mathematical  analysis  of  the  longitudinal  impact  of  a_  round 
head  short  bar  on  an  infinitely  long  bar. 

This  discussion  follows  that  of  W.  A.  Prowse  [1]  (see  Figure  12). 
Let  the  round-head  anvil  strike  the  plane  surface  with  velocity 
v.  Let  P (t)  be  the  total  indentation  of  the  anvil  and  surface  com- 
bined. Then,,  for  steady  pressure  p0  the  Hertz  law  [13]  gives 

3/2 


p = k2g 


(3.57) 
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where 


2E 

3(1  + v ) u - v ] 


(3.58) 


E and  v are  the  Young's  modulus  and  Poisson  ratio  of  the  anvil 
and  surface  while  r is  the  radius  of  the  anvil  head. 

Prowse  considers  a section  of  the  anvil  of  length  Ax, 
traversed  by  the  wave  with  velocity  c in  time  At.  The  resul- 
tant force  acting  on  the  face  AB  of  this  section  is  equal  to 
the  time  rate  of  change  of  momentum  of  the  section. 

The  resultant  force  is 


(3.59) 


while  Prowse  considers  the  time  rate  of  change  of  momentum  as 


HF  (ApA4  HF}  “4  pcAAt^l 


.d^p 


(3.60) 


where  A is  the  cross-sectional  area  of  the  anvil. 

Thus,  equaling  (3.59)  and  (3.60) 

"HF  (P  * \ P cA  h|)  * 0 

Using  (3.57)  this  yields 

k2P3/2  * -jpc  Ag|  = cx  (3.61) 

We  may  use  either  initial  or  final  conditions  for  the 
determination  of  the  constant  of  integration.  First  we  follow 


52 


Prowse  and  designate  the  maximum  indentation  by  Thus  when 


(3.62) 


P = P- 


and  it  follows  from  (3. 61) that 


k2(Eli3/2  * p3/2) 


(3.63) 


Since  the  initial  indentation  is  zero,  the  solution  is 


pcA  r W1  

IkloP2  ° 1 - 


(3.64) 


Prowse  left  the  result  in  this  form  and  evaluated  the  integral 
numerically.  For  our  purposes  it  is  better  to  evaluate  the 
integral  analytically  before  starting  numerical  computations. 


Define 


s - KpTp 


(3.65) 


where  v is  the  initial  velocity.  Then  we  have  from  (3.64) 

* -#  ■ tan'x  (¥i)  (3-66) 
From  the  Prowse  condition 


7^ 


= v for  3=0 
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we  obtain  an  expression  for  the  limiting  indentation  in  terms 

of  the  initial  velocity,  from  (3.57) 


0 


1 


3/2 


(3.67) 


Note  that  the  limiting  stress  is 


a -^1 

al  K 


3/2 


Consequently,  by  (3.67)  and  (3.68) 


(3.68) 


= Pcv  (3.69) 

which  is  the  stress-velocity  relation  of  the  theory  of  impact  due 
to  St.  Venant  [12],  [13]. 

It  is  possible  to  construct  stress-time  curves  based  on 
equation  (3.66)  and  Table  2.  This  has  been  done  in  Figures  13 
through  28  for  steel  hammers  with  hemispherical  ends.  The  critical 
values  for  these  figures  are  given  in  Table  3. 

L. ..ncothe  hammer  length 

h. .... .the  height  of  drop 

2L 

t^o.o.obase  time  = (time-length  of  corresponding 

St.  Venant  pulse) 

To . . o o o time-length  of  pulse 
cT.oo oo liner  pulse- length 
a o o o . omaximum  stress  attained 
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The  velocity  of  the  hammer  at  time  of  contact  and  the  theoretical 
limiting  stress-,  are  independent  of  the  length  of  the  hammer  and 
are  given  as  a function  of  the  height  of  drop  in  Table  3. 

For  these  computations  the  following  values  were  assumed: 

E * 2„073  x 1012  dynes/cm2 

v = 3 x 10‘ 1 

3 

p = 7„7  gm/cm 

Hammer  diameter  - 2„5  cm 


Figure  12 

Round-Head  Hammer  Impact 


TABLE  3 

Critical  Values  for  Figures  13  through  28 


Fig. 

L 

(cm) 

h 

(cm) 

, >1  . 

(micro) 

sec) 

T 

(micro 

sec) 

cT 

(cm) 

°n  2 
(KgF/ cin  ) 

10 

25 

5 

96.5 

140 

72.5 

348 

11 

25 

10 

96.5 

136 

70.5 

514 

12 

25 

15 

96.5 

134 

69.2 

644 

13 

25 

25 

96.5 

130 

67.5 

850 

14 

50 

5 

193 

239 

124 

405 

15 

50 

10 

193 

233 

121 

575 

16 

50 

15 

193 

231 

120 

708 

17 

50 

25 

193 

228 

118 

914 

18 

75 

5 

290 

336 

174 

409 

19 

75 

10 

290 

331 

171 

578 

20 

75 

15 

290 

328 

170 

708 

21 

75 

25 

290 

325 

168 

914 

22 

100 

5 

386 

432 

224 

409 

23 

100 

10 

386 

427 

221 

578 

24 

100 

15 

386 

424 

220 

708 

25 

100 

25 

385 

421 

218 

914 

1.0 
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Stress-Time  Relationships 


Stress-Time  Relationships 
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Stress-Time  Relationships 
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Stress-Time  Relationships 


~ 64 

j 


1„0 
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Stress-Time  Relationships 
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Stress-Time  Relationships 


70 


Stress-Time  Relationships 


1.0 
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Str 
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Stress-Time  Relationships 
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SAMPLE  COMPUTATION  AND  DISCUSSION  OF  TEST  RESULTS 
FOR  NATURAL  RUBBER  (HEVEA) 


SECTION  IV 

Sample  Computation  and  Discussion  of  Test  Results 
for  Natural  Rubber  (Hevea) 


We  consider  the  case  of  normal  (approximate)  velocity  of 
50  cm/sec.  The  material  is  natural  rubber. 

The  raw  data  obtained  by  averaging  readings  from  film  strips 


for  ten  separate  runs  are  as  follows 

Calibration:  (initial  length)  185.3 


Average  drum  speed:  121,17  rpm 


Readings : 


~3 

Lengthwise  (.10  in.) 

Transverse 

(Calibration) - (Width) 

0 

185.3 

0 

100 

179.3 

6.0 

200 

172.9 

12.4 

300 

167.0 

18.3 

400 

160.9 

24.4 

500 

154.8 

30.5 

600 

149.6 

35.7 

700 

144.1 

41.2 

800 

139.2 

46.1 

900 

134.9 

50.4 

1000 

131.2 

54.1 

1100 

127.8 

57 . 5 

1200 

125.3 

60.0 

1300 

123.4 

61 .9 

1400 

122.4 

62 .9 

1500 

123.1 

62.2 

1600 

124.0 

61.3 

1700 

125.7 

59.6 

1800 

128.9 

56 .4 

1900 

132.7 

52,6 

2000 

137.0 

48  • 3 

2100 

142.0 

43.3 

2200 

2300 

148.0 

154.0 

37.3 

31.3 

2400 

2500 

2600 

160.7 

167.9 

175.2 

24 . 6 
17.4 
10.1 

2700 

2745 

2800 

182.6 

185.3 

190.0 

2.7 

0 

-4.7 
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These  results  are  then  plotted  and  a smooth  curve  drawn  through 
the  points  (Figure  29,  30,  31), 

The  scale  of  the  curve  is  changed  into  the  physically— 
meaningful  parameters,  time  and  displacement,  through  the  following 
computations. 

The  inside  diameter  of  the  drum  is  15.06  inches,  its  average 
rotational  speed  is  121.17  rpm.  Hence  the  time  corresponding  to 
a given  horizontal  distance  h in  inches  on  the  film  is 

60  h 

1 TsmmsmBff 

hence  the  total  time  of  contact  is 

T~  HI7I7  x*T5;T)'Sy  * 2.87244  x 10  2sec 

The  length  of  the  specimen  is  1.27  cm.  The  calibration  measure- 
ment yields  185.3.  Thus  our  change  of  scale  on  the  displacement 
is  given  by 

q « 1.27  . 

S 185.3  d 

We  divide  the  abcissa  of  the  graph  into  10  divisions  of 
length  2.87244  x 10-3sec  and  take  the  values  of 

5(k  ) * S($)  » Sk 


from  the  graph.  We  have 


Displacement  (cm) 
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S (cm) 


Figure  29 

Displacement "Time  Curves 

_3 

Nominal  velocity  25  cm/secj  co  = 3.088x10  sec 
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S (cm) 


Figure  30 

Displacement-Time  Curves 

-3 

Nominal  velocity  50  cm/ sec;  co  = 2.872x10  sec 


Displacement  (cm) 
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Time  (sec) 
Figure  31 


Displacement-Time  Curves 

-3 

Nominal  velocity  75  cm/sec;  co  * 2.645x10  sec 
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k 


s 


k 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0,0 

0.1330 

0.2563 

0.3530 

0.4147 

0.4297 

0.3934 

0.3228 

0.2262 

0.1179 

0.0 


We  wish  to  fit  this  data  to  a function  of  the  form 


s(t)  = C[e"nt  + | e“(p  + iq)t  +|e‘(p  “ iq)t] 


where 

,2.2 

k - -1  ♦ i -is.  . g - P ) 

2qp 

This  form  was  chosen  to  insure  that  the  displacement  and  accelera 
tion  are  zero  fot  t = 0. 

We  use  the  procedure  outlined  by  Whittaker  (see  the  quotation 
following  this  discussion),  with  the  ordinary  method  of  least 
squares  to  find  the  values  of  the  constants  Dg,  D-^,  Dq  which 
most  nearly  satisfy  the  set  of  equations 

s3  + D2s2  + ^lsl  + D0S0  “ 0 
s4  * D2s3  + D1s2  + D0S1  * 0 
s10  + D2S9  * D1s8  + ®0S7  “ 0 
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This  procedure  leads  to  the  system  of  equations 

7.338 ,1605  D2  + 8.181,983,5  Dx  + 8.235,692,7  DQ  » -5.816,025,7 
8.448,673,3  T>2  + 8.747,357,1  -+  8.181,983,5  Dx  * -7.338,160,5 

8.709,471,2  T0  + 8.448,673,3  - 7.338,160,5  Dx  * -8.107,794,3 

which  has  the  solution 

D2  - - 2.068,220,275 

Dx  * 1.350,660,495 

,D0  = -0.205,  225,635,7 

We  now  solve  the  equation 

x3  + D2x2  + Dxx  + Dq  - 0 


for  its  roots 


*1  * e 


-no) 


- (p  + i<aW 


. xo 

* , 

and  find,  after  obvious  calculations,  -that 


- (p  - iq)o> 


n » 5.337,880  x 102,  p - 8i,  771,98,  q - 1.107,459  x 102 
Now  that  we  have  n,  p,  and  q,  we  caiu  "write 

s(t)  * C[e'nt  - e'ptcos  qt  * (n  * ^qp  P~  0 Ptsin  qt] 


and  determine  the  value  of  C which  b-est  satisfies  the  equations 
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Mk 


s,  * G[e‘^  - e^Gos  q^ 


,2^2  2v 

+ (n  j_a_ i-EJ  „*p 


iqp 


e ‘ sin  qojJ 


s,  - CCe'2,w  - e-2P“ooS  2q{0  + (,P2  g?>.  e'^sin  2<W] 


on>saaaae«Daac«a»ao 


ooBeooocoooaoo 


’10 


G[e 


■io™  - -*10':>“cosloqu+  <n.2 


■ ■aaifacootato 


2qp 


lOpw 


sin  lOqw] 


We  again  use  the  method  of  least  squares,  finding  this  time  that 


C * 0„002, 915, 748 


The  hereditary  constants  a,  Hq  and  H are  obtained  directly  with 

the  use  of  the  secant  modulus  for  this  material  and  maximum  de- 

7 2 

formation  E = 2„81  x 10  dynes /cm  . Hence  we  have 

la 

a = n + 2p  = 5 „ 513 ,320  x 102 

H = ^2.  (p2  -*•  q2)n  - E2=  lc  626, 555  x 107  dynes/cm2 
o cm0 

H - (2pn  + p2  + q2)  - Hq  - E2  = 3,662,901 

7 2 

x 10  dynes/ cm 


We  have  two  checks  on  the  accuracy  of  our  computations.  The 
first  is  a comparison  of  the  initial  velocity  obtained  from  the 
graph;  vQ  = 47,5  cm/sec  with  the  value  obtained  from  the  corn- 

put  at  ion 


dsn 

HP  tO 


C[p-n 


n 


+ q 

“2p 


2 ■ 


-E— ] = 47„849  cm/sec 
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The  second  check  is  a straight  comparison  of  the  experimental 
values  s^  with  the  values  s(kto)  obtained  from  the  equation,, 

The  latter  check  yields  an  error  of  less  than  4 percent  over  the 
entire  range  of  time.  Further,  all  values  were  rounded  off  only 
after  the  computation  was  completed  rather  than  at  any  inter- 
mediate stages. 

Whittaker's  summary  (Reference  8)  of  the  Prony  method  of 
approximating  a function  with  a finite  number  of  exponentials  is 
quoted  below.  The  notation  has  been  altered  to  agree  with  that 
of  this  report  as  much  as  possible. 

"Although  Prony's  method  is  more  than  a century  old, 
it  does  not  appear  to  be  widely  known  or  to  have  found 
its  way  into  any  text-book;  and,  as  his  original  paper 
is  perhaps  not  accessible  to  many  readers,  I may  be 
justified  in  giving  here  a brief  notice  of  it." 

"Suppose  that  s(t)  is  given  numerically  for  a cer- 
tain range  of  values  of  t.  Take  any  set  of  values  of 
t equally  spaced  within  this  range,  say  t - 0,  o>  , 2a)  , 

3(0  , 4 o),  ...,  and  let  the  corresponding  values  of  s(t) 
be  sQ,  s1#  s2,  s3,  ....  Now  if  s(t)  could  be  represented 
exactly  in  the  form  of  a sum  of  k exponentials,  say, 

Pept  + Qeqt  + Rert  + . . . + Vev"t, 
then  s ( t ) would  satisfy  a linear  difference-equation  of 


the  form 
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AVfc  + BVk-l  + CVt-2  + •••■  + Msn  " 0 

where  the  roots  of  the  algebraic  equation 

Atk  + Bi^1  + Ctk“2  + + M = 0 

would  be 


"Prony's  method,  which  is  based  on  this  fact,  is 
to  write  down  a set  of  linear  equations, 

Ask  + Bsj.^  ♦ Csk.2  + + Ms0  - 0 

As^.+  ^ + Bs^  + Gsk-^  + .....  + Ms^  = 0 

Ask+2  + Bsk+i  + Csk  + + MS2  “ 0 

Ask+3  + Bsk+2  + Csk+1  + + Ms3  " 0 

where  the  quantities  Sq,  s^,  Sg,  Sj,  are  known, 

since  s(t)  is  a known  tabulated  function,  and  by  the 
ordinary  method  of  Least  Squares  to  find  the  values 
of  A,  B,  C,  ...,  M,  which  most  nearly  satisfy  the 
equations;  then  with  these  values  of  A,  B,  C,  ...,  M, 
to  form  the  algebraic  equation 

Ask  + Bsk~X  + Csk~2  + ...  + M - 0, 

and  find  its  roots;  these  roots  will  be  ep  , eq  , er  , 
...ev  , and  thus  p,  q,  r,  ...,  v,  are  determined. 
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Knowing  p,  q,  r,  , . . „ v,  we  have  a set  of  linear  equations 
to  determine  the  coefficients  P,  Q,  R,  . . „ , V,  and  these 
also  are  to  be  solved  by  the  method  of  Least  Squares,  " 


Discussion  of  Test  Results  for  Natural  Rubber  (Hevea) 

The  data  obtained  from  a series  of  tests  on  a natural  rubber 
27-30  durometer  hardness  were  analyzed,  by  use  of  the  procedure 
outlined  above,  for  three  initial  speeds  of  the  striking  bar. 

It  was  found,  as  might  be  expected,  that  the  phenomena  were 
essentially  viscous  for  nominal  speeds  of  approximately  25,  50, 
and  75  centimeters  per  second.  Indeed,  in  the  case  of  the  highest 
initial  speed  (and,  correspondingly,  the  greatest  distortion), 
the  hereditary  effect  was  so  small  as  to  be  obscured  by  inaccuracies 
of  the  data. 

The  pertinent  results  of  the  investigation  are  given  in 
Tables  4 and  5,  which  follow.  For  the  analysis  which  was  made  under 
the  assumption  of  a viscous  stress-strain  laws 


0 - De  * at 

s (t)  - Ce  ptsinqt 


(4.1) 


We  list  in  Table  4 the  parameters,  p,  q,  C,  D,  and  y;  as  well  as 
the  time  of  contact  T,  the  actual  initial  velocity  v f and  the 
maximum  strain,  emo 

The  slight  hereditary  effect  displayed  by  the  hevea  indicated 
that  it  would  be  pointless  to  carry  the  analysis  further  than 
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heredity  of  the  first  degree  with  residual;  which  is  characterized 
by 

a “ £e  + / t [H0  + He“a(t''T,  ||  dx  (4.2) 

o 

The  parameters  E,  Hq,  H,  a,  C,  n,  p#  q,  of  (4.2),  are  listed  in 
Table  5.  Note  that  E is  the  static  secant  modulus  for  natural 
rubber. 

We  finally  remark  upon  an  interesting  connection  between 
these  two  analyses.  Comparison  of  Tables  4 and  5 reveals  that, 
for  a given  initial  velocity,  we  have  approximately 


D = E + Hc 


Thus,  it  can  be  concluded  that 


f He 


-a(t-x) 


de 

cF 


dx 


for  the  hevea  specimen. 
Since  we  know  that 


lim 
a— > oo 


* -a(t-x)  de 


3F 


uT 


de 

cF 


we  expect  that  the  approximate  relation 
H * aa 


will  hold. 
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TABLE  4 

Results  - Viscous  Case 


Nominal 

Velocity 

cm/sec 

sec 

vo 

cm/ sec 

e 

m 

P 

q 

C 

cm 

D 2 
dyne  s / cm 

, Y 
dyne- sec 

cm2 

25 

3. 09x10" 2 

23.7 

0.172 

5.15 

102 

0.233 

3.86xl07 

4.02xl04 

50 

2. 8 7x10" 2 

47.8 

0.338 

5.05 

110 

0.453 

4 . 47xl07 

3. 75xl04 

75 

2. 65xlO"2 

77.3 

0.481 

6.08 

119 

0.650 

5. 27xl07 

4. 52xl04 

TABLE  5 

Results  ° Hereditary  Case 


Nominal 


H 


H 


a 


Velocity 
cm/ sec 

n 

P 

q 

E £ 
dynes /cm 

o 2 

dynes /cm 

dynes /cm2 

“1 

sec  cm 

25 

642 

5.36 

102 

2.44xl07 

1. 36xl07 

2. 62xl07 

653  6.02x10"^ 

50 

543 

8.77 

111 

2. 81xl07 

1. 63xl07 

3, 62xl07 

551  2. 9 2x10" ^ 

These  results  are  shown  graphically  in  Figs.  32  through  39 


8x10 
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True  Strain 


o)2(i)3a)4a)5a)6(o  7(0  8a>9o)  IGo 

Time  (sec) 

Figure  33 

Hereditary  Function 

-3 

Nominal  velocity  25  cm/ sec;  co  * 3.088  x 10  sec 


0) 


2o)  3co  4ti)  5(0  6o)  7(0  8o)  9o)  10to 

Time  (sec) 

Figure  34 

Hereditary  Function 

Nominal  velocity  50  cm/secj  co  = 2.872x10  v sec 


Stress  (dynes/cm  ) 
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Stress  versus  Strain  - Viscous  Analysis 
Nominal  velocity  25  cm/ aec 


Stress  (dynes / cm 
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(N 


a 


Figure  36 


Stress  versus  Strain  - Viscous  Analysis 
Nominal  velocity  50  cm/ sec 


Stress  versus  Strain  - Viscous  Analysis 


Nominal  Velocity  75  cm/ sec 


Strain 


Figure  39 


Stress  versus  Strain  - Hereditary  Analysis 
Nominal  velocity  50  cm/sec 
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At  present;  a series  of  tests  has  been  completed  on  the 
Butyl-25  and  Neoprene  GNA.  The  computation  of  the  results  is 
now  underway.  These  tests  will  be  followed  by  a aeries  on 
1.  Nylon 
2«  Polythene 
3,  GRS  artic. 

The  results  of  these  tests  will  appear  in  a final  report  in 
June  1953* 
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APPENDIX  A 


MATHEMATICAL  ANALYSIS  OF  THE  LONGITUDINAL  IMPACT 
OF  A BALL  ON  AN  INFINITELY  LONG  BAR 


APPENDIX  A 

Mathematical  Analysis  of  the  Longitudinal  Impact 
of  a Ball  on  an  Infinitely  Long  Bar 


A steel  bar  of  mass  mf  radius  r,  and  velocity  vQ  strikes 
the  end  of  an  infinitely  long  bar.  We  wish  to  find  the  pressure 
at  the  end  of  the  bar  as  a function  of  time. 

We  take  the  stresses  positive  in  compression  and  the  dis- 
placement positive  directed  into  the  bar. 

For  an  arbitrary  force  F(t)  applied  uniformly  over  the 
end  of  the  bar: 

a(x,  t)  * F(t  - £) 

where 


Since  u * u(t  * -)  we  have 

/ r-du  _ n du 

a(x,  t)  = Pc-SF 

The  complete  displacement  due  to  the  impact  is 
d - P + u 

where  fi  is  the  total  indentation  and  u is  the  elastic  deformation 
due  to  the  force  F(t).  The  velocity  is 
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By  the  law  of  impulse  and  momentum,  we  have  at  x * 0 

! t 

T - - £ / F(t)  dt 

o 

thus 


yc  - |/F(t)dt 


Assume  the  Hertz  law  [13],  i.e.. 


F - 


3/2 


2E 


k2  ‘ 3(1  + v) (1-  v) 


rl/2 


thus 


(p3/2)  + p3^2  - o 


This  differential  equation  is  not  easily  integrated.  Numerical 
solutions  for  boundary  conditions  of  the  type 


P • 0] 

l for  t » 0 

v " voj 

are  available,  however,  for  specific  values  of  vQ  and  the  con- 
stants of  the  equation. 

Assume 

k - 1.73  x 1012  d>rn|jK  m - 365.1  gms 

2 cm' ' 

pc  - 4.05  x 106  321S®£  A - 5.07  cm2 

cm 
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then  for 


where 

b = —4;  d = ™ 

pcA*  m 

we  have 

b * 8.44  x 104  (cm'^sec)"1 


d * 4.75  x 109  (cm1/2sec2)~1 


Initial  conditions: 


t * 0S  3 - Oj  * vQ 


For  m = 365.1  gms(  five  cases  were  considered: 


TABLE  VI 

Ball  Impact  Values 


vo 

(cm/ sec) 

T 

(|i  sec) 

3 

*max 

(mm) 

F 

tmax 

(x  10adynes) 

a 

max  c) 
(Kg  F/ cnr ) 

45 

245 

0.0224 

5.80 

1166 

90 

19  2 

0.0367 

1.22 

2453 

135 

176 

0.0484 

1.84 

3704 

180 

164 

0,0598 

2.53 

5088 

220 

157 

0.0699 

3.20 

6435 

The  graphical  solutions  are  given  in  Figures  40  through  44 „ 


(See  Table  6 ) 


Ball  impact  - v =90  cm/sec 
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APPENDIX  B 


Rubber  Specimen  Recipes 


The  1/ 2-inch  by  1/2-inch  rubber  cylinders  were  fabricated  in 
the  Rubber  Laboratories,  Armour  Research  Foundation  of  Illinois 
Institute  of  Technology#  They  were  made  in  an  eight  cavity  mold 
of  the  following  rubber  compounds: 

1.  Hevea 

2.  Hycar  OR- 15 

3.  Neoprene  GNA 

4.  GRS  Artie  (85  Butadiene/15  Styrene) 

5.  Butyl- 2 5 

The  compounds  were  prepared  on  a two-roll  6 x 10-inch 
laboratory  rubber  mill  with  a front  roll  speed  of  30  fpm  and  a 
speed  ratio  of  back  roll  to  front  roll  of  1„4  to  1,  Cold  water 
was  circulated  through  both  rolls  during  the  mixing  of  the  GRSff 
Butyl,  Neoprene,  and  Hycar  0R-15  compounds „ No  water  was  cir- 
culated during  the  mixing  of  the  Hevea. 

The  formulas  for  the  compounds  are  given  in  Table  7 and  a 
summary  of  their  physical  properties  is  given  in  Table  8. 
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TABLE  7 

RUBBER  SPECIMEN  FORMULAS 


GRS  (Artie) 


GRS  (artic)  100 
Reogen  10 
p-33  25 
Zinc  Oxide  XX78  5 
Stearic  Acid  2 
Agerite  Stalite  1 
H Teads  3 
D.PoG.  1 


Neoprene 


Neoprene  GNA 

100 

L„C.  MgO 

3 

P-33 

25 

Circo  Lt.  Process  Oil 

15 

Neozone  D 

1 

Stearic  Acid 

1 

Zinc  Oxide  XX78 

5 

Hycar  OR-15 


Hycar  OR-15  100 
p-33  25 
Dioctyl  Adipate  30 
Zinc  Oxide  XX78  5 
Stearic  Acid  1 
Agerite  Stalite  1 
H Teads  3 
D.P.G.  1 


Butyl  25 


Butyl  25  100 

P-33  25 

Zinc  Oxide  XX78  5 

Mineral  Oil  (White)  1 

Lime  (CaOH)  5 

Sulfur  2 

Thiomex  1 

Polyac  0.5 

Akroflex  C 2 


Natural  Rubber 


Smoked  Sheets 

100 

Reogen 

10 

Zinc  Oxide 

5 

Stearic  Acid 

2 

Agerite  Stalite 

1 

P-33 

25 

Methyl  Teads 

3 

D.P.G. 

1 
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TABLE  8 

SUMMARY  OF  PHYSICAL  PROPERTIES  OF 
RUBBER  COMPOUNDS  GIVEN  IN  TABLE  7 


Stock 

320°F  M500 

T 

E 

S Shore  "A"  Hardness 

Hevea 

10 

600 

2300 

850 

15 

38 

20 

710 

2880 

710 

10 

38 

30 

725 

2550 

670 

8 

38 

GRS  (Artie) 

30 

215 

270 

565 

10 

38 

Hycar  OR- 15 

30 

195 

1280 

945 

13 

31 

Neoprene 

30 

965 

2620 

840 

15 

44 

Butyl  25 

30 

365 

285 

0 

43 

M500  - 

Modulus 

at  500  per 

cent  Elongation 

T - 

Tensile 

Strength  (psi) 

E - 

Elongation  (per  ce 

nt) 

S - Permanent  Set  at  Break 
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